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We investigate the entanglement entropy between two subsets of particles in the ground 
state of the Calogero-Sutherland model. By using the duality relations of the Jack symmet- 
ric polynomials, we obtain exact expressions for both the reduced density matrix and the 
entanglement entropy in the limit of an infinite number of particles traced out. From these 
results, we obtain an upper bound value of the entanglement entropy. This upper bound has 
a clear interpretation in terms of fractional exclusion statistics. 

PACS numbers: 



I. INTRODUCTION 

Entanglement properties of quantum many-body systems have been attracting much attention 
in quantum information theory and condensed matter physics. The entanglement entropy (EE), i.e., 
the von Neumann entropy of the reduced density matrix of a subsystem, is a measure to quantify 
how much entangled a many-body ground state is. Recently the EE has been used to investigate 
the nature of quantum ground states such as the quantum phase transition and topological order 
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Sj. When we study the entanglement properties in many-body systems, exactly solvable 



models in one dimension such as the harmonic chain [6|] , the XY spin chain in a transverse magnetic 
field l|, Q, and the Affleck-Kennedy-Lieb-Tasaki model 9, 10, [ill, [3l serve as a laboratory to 
test the validity of this new concept. The relation between the EE in solvable models and the 
conformal or massive integrable field theories is extensively discussed in Refs. 

In this article we study the EE of the ground state of the Calogero-Sutherland (CS) model 

nn 

[151. Il6l|. The CS model is a quantum integrable model with inverse-square interactions on a circle. 
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An infinite number of conserved quantities which characterize the integrable structure of this model 
have been constructed in a Lax form 17]]. Although it is usually a formidable task to compute 
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8K one can derive exact expressions for 



2l[. This is an important feature of this 



the correlation functions even in the integrable mode, 
the dynamical correlation functions in this model 1^, 
model which distinguishes itself from the other integrable models. Another interesting aspect of this 
model is a connection with the fractional statistics in low dimensions. In fractional quantum Hall 

n 

systems, the ground state wave function is given by the Laughlin state [231], and its excitations have 
fractional charges. Similarly, the ground state of the CS model is described by the Jastrow-type 
wave function and its excitations are also quasiholes with fractional charges. Then we can identify 
the CS model as a canonical model to study the exotic properties of the fractional statistics in low 



dimensions. It s 
studied recently 



hou 
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d be noted here that the EE of the Laughlin state itself is also extensively 
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26|. 



We consider the EE between two subsystems in the ground state of the CS model. Let us first 
explain how to partition our total system into two subsystems. There are mainly two possible 
ways to partition the system under consideration. One way is to divide the system into two spatial 
blocks the other to divide the A^-particle system into an L-particle block and an (A^ — L)-particle 
block. They are called a spatial partitioning and a particle partitioning, respectively. In this 
article, we focus on the latter. As the EE between two spatial regions in the fractional quantum 
Hall states can extract a topological quantity such as the total quantum dimension [24 1, the EE 
based on the particle partitioning in the CS model reveals a new aspect of low-dimensional systems 
with fractional exclusion statistics. First we consider the L-particle reduced density matrix of 
our system. By using duality relations of the Jack polynomials, we can formally obtain the exact 
expression for the reduced density matrix. Although we have the exact form of the reduced density 
matrix, it is difficult to evaluate the eigenvalues since there are many off-diagonal elements. Then 
we consider the thermodynamic limit and find that a great simplification occurs in this limit. 
We should note here that what we mean by the thermodynamic limit is {N — L) ^ 00 limit, 
where (A^ — L) is the number of particles traced out. It is slightly different from the usual sense 
such as A^ — > 00 with fixed L/N. Finally, we focus on the upper bound value of the EE. In the 
thermodynamic limit, we can approximate the reduced density matrix by a maximally entangled 
state and hence we can evaluate the upper bound by counting the allowed Young tableaux in the 
duality relation. The upper bound value is estimated as S^Y^'^ = log (^^(^~^^)+-^^ and has a clear 
interpretation in terms of exclusion statistics [^l- We also find that the subleading term of the 
EE is independent of the total number of particles N. 
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The organization of this article is as fohows. In Section |TT1 we wih introduce some basic concepts 
in the CS model used in later sections. Section Hill is the main part of this article. We will calculate 
the reduced density matrix in the CS model and show that it becomes very simple if we take a 
thermodynamic limit. Then we will be able to obtain the EE in this limit and to estimate the upper 
bound of this EE. We will discuss the physical interpretation of this upper bound. Section HVl will be 
devoted to summary and discussions. In AppendixEl we will analyze the EE in the thermodynamic 
limit more in detail than Section IIIII 



II. CALOGERO-SUTHERLAND MODEL AND JACK SYMMETRIC POLYNOMIALS 

A. Calogero-Sutherland model 

We introduce a precise definition of the CS model. The CS model describes the interaction of 
N particles on a circle of length / and the Hamiltonian is given by 

^tt2 9x|+^sin^(f(x,-x,))' 

where xj (0 < xj < I) are the coordinates. Here it is convenient to introduce new coordinates on a 
unit circle zj = exp(^Xj). Using these new variables, the exact ground state of Hcs is given by 
the Jastrow-type wave function as 

1 ^ 

All the excited states of this model can also be obtained by multiplying certain symmetric poly- 
nomials to ipQ as 

■0a(^1,22, ■■■,Zn) = Px{zi,Z2, ...,ZN;P)'lpo{zi.Z2, .■.,Zn). (3) 

The symmetric polynomial in Eq. ([3]) are called the Jack symmetric polynomials and characterized 
by partitions A. The partition A is a sequence A = (Ai, A2, Ar, ...) of non-negative integers in 
decreasing order: Ai > Ao > ... > Ar > ... . Let us introduce some terminology. We use the 



notation of Macdonald 29|. Every partition has a corresponding Young tableau which graphically 
represents a partition (see Fig. The non-zero Aj are called the parts of A. The number of parts 
is the length of A, denoted by Z(A) and the sum of the parts is the weight of A denoted by |A| and 
explicitly written as |A| = Yl\=i The excitation energy is also characterized by the partition as 

N 

^a = ^(t)'E^'(^)' (4) 
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(a) ■■ -0 (1 0) (1 0) (1 0) (1 0) (1 0) 

(b) ■■ -0 9 (1 0) (1 0) (1 0) (1 0) ■ 

3 quasiholes ; 



FIG. 1: (a) The ground state configuration of tlie quasi-momenta a,t (3 — 3 witli 5 particles, (b) Tlie excited 
state with 4 particles and 3 quasiholes obtained by removing one pariticle from the Fermi Sea. 

where the quasi-momentum ki{X) = Aj + — i)- The set of quasi-momenta is subject to the 

exclusion constraint ki — /cj+i > (3. In the ground state, the configuration of the quasi-momenta 
is given by A;j(0) = P{^^y^ — i) and this configuration is schematically shown in Fig. [TJ(a). We 
call this configuration the Fermi sea. In Fig. [H a particle can be identified by one 1 followed by 
P — 1 zeros and a quasihole by one 0. Therefore, if we remove n particles from the Fermi sea, /?n 
quasiholes are created in the Fermi sea (see Fig. [11(b)). We should note here that the coupling /3 
has been assumed to be a positive integer for the sake of simplicity in this article. However, in 
principle, we can extend this correspondence at any positive rational coupling (3 = p/q 
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B. Jack symmetric polynomials 

Let us turn to focus on the mathematical aspects of the Jack symmetric polynomials. The Jack 
symmetric polynomials are mutually orthogonal with respect to the following scalar product on 
the ring of symmetric polynomials in N indeterminates zi,...,zm'- 

dzi f dzN 



{f^9)'N= ^ J ^/(^l,^2,-",^;Af)5'(^l,^2,"-,^Ar)|V'o(2l,^2,-",^;Af)| • (5) 



The normalization of the ground state wave function ipQ is defined as J\f{f3,N) = (1, 1)^ and its 



explicit form is given by j^TNkn- The explicit orthogonality relation for the Jack polynomials is 



given by 

where s = is a box on a Young tableau identified by its coordinates 1 < i < 1{X) and 

1 < j < Aj. The notations a{s), l{s), a'{s) and l'{s) are summarized in Fig. [2j It is well known that 
classical families of symmetric polynomials can be obtained by specializing the coupling /? of the 
Jack symmetric polynomials. For /? = 0, 1, 2, and oo, the Jack symmetric polynomials are reduced 
to the monomial symmetric, the Schur, the zonal, and the elementary symmetric polynomials, 



respectively 



29|. 
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FIG. 2: The young tableau corresponding to the partition A. The box in the tableau indicates s = ihj)- 
a{s),l{s), a' (s) and l'{s) are arm-length, leg-length, arm-colength and leg-colength, respectively. 

III. REDUCED DENSITY MATRIX AND ENTANGLEMENT ENTROPY 



In this section, we consider the reduced density matrix and the entanglement entropy for any 
subset of L particles in a system of N particles in the state The L-particle reduced density 
matrix, being normalized, is defined as 



p{wi, ...,wl;zi, ■■■,zl) 

1 f dzL+1 



dzjy 



Here the partial trace is taken over the variables z^+i, ztv- To calculate the EE, it is use- 
ful to introduce a trace in a complex integral form. The trace of any L-particle operator 
A(Wi, ...,WL', zi, zl) is defined by 

Since the reduced density matrix ([7]) is normalized, Tr [p] = 1. Similarly, the trace of the product 
of any L-particle operators A{wi, ■.■,wl', zi, zl) and B{wi, ■.■,wl', zi, zl) is defined by 

dwL f dzi f dzL 



J 2Txiwi J 2'KiwL J 27rizi / 



„ . r ^ . r ^ . -A{wi,...,wl;zi,...,zl)B{zi,...,zl;wi,...,wl), 

and the EE is defined by Sn^l = — Tr [plogp]. To obtain the explicit form of the reduced den- 
sity matrix, it is convenient to rewrite Eq. ([7]) by using the ground state wave functions of the 
subsystems, iIjq{zi, Z2, zl) and V'o(-ZL+i, -^Af), as 



p{wi,...,wl]Zi,...,zl) 



1 L\{N-L)\ 
M{(3,N) Wl 



n 

i=l 



WiZi 



ipoiwi, ...,wl)iPo{zi, Zl) 



dzL+i 



L N 



f^^Il n a-^^Z,)''{l-mZ,)('\MzL+l,•■.,ZN)\' 



2-nizL+i J 2mzN V r , i 

1=1 ]=L+1 



(9) 
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FIG. 3: Young tableaux within the shaded region are aUowed in the expansion formula pT|) . 



Recalling the definition of the scalar product (l5|), Eq. ([9]) can be rewritten again as 



I I L N 

^^!Q{wi,...,WL)^Q{zi,...,ZL)i]\ n (1 



N 



(10) 



where ^o{zi^ ...^ zl) = (Hi^i -^«)~^^^~'^^^^V'o(-2i) -^l)- The next thing to do is to compute the 
scalar product in Eg. (fTO|) . Let us now introduce the following duality relation to carry out our 
calculation 
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WWO- + XiVj) = ^P\{xi,x2,...,XN;l3)Px'{yi,y2,--,yM;l/P)- 

i=lj=l A 



(11) 



Here, the conjugate partition A' is a transpose of the Young tableau A and partitions A are summed 
over the Young tableaux which satisfy Z(A) < N and /(A') < M (see Fig. ([3])). The duality relation 

Eq. (jlip plays a crucial role to simplify the reduced density matrix (jlOp . We shall explain the 

(k) 

procedure of the calculation in more details. First, we introduce dummy variables Zj,[L + l< 
j < N, 1 < k < (3). Secondly, we expand ni^in(jfc)(l ~ using the duality relation 



(jlip . Here, (j, k) runs from (L + 1, 1) to {N,P). Finally, we set the dummy variables z 
(1 < < /?). We can summarize the above as the following expansion formula: 



(fe) 



N 



/3 



/3 



n n (1 



^ Px{zi, ZL;P)Py{-ZL+l,..; -ZL+1, -ZN, -zn] 1//?), (12) 



i=li=L+l A 

where partitions A are summed over those that satisfy 1{X) < L and /(A') < P{N — L). Here we 
have also assumed that the coupling /3 is a positive integer. The above formula has a clear physical 
interpretation as a superposition of the intermediate states consist of L particles and [3{N — L) 
quasiholes. 
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Next, we try to rewrite Py with coupling in Eq. (jl2p in terms of Px with (3. It is also 
well known that the Jack symmetric polynomials can be expressed as polynomials in power sums 
Pn = zf-. We give as examples the expressions up to |A| = 3: 



^(1) =Pi, 

P{2) = J^P2 + J^pI 

P{i,i) = -hp2 + \pI, 



(3) = (l+/3)(2+/3)^'3 + (l+/3)(2+/3)^'2Pl + (i+/3)(2+/3)^'l' 
^^(2,1) = -TW^3 + ^P2Pi + t^pI (13) 

, ^(1,1,1) = \P3 - \P2PI + \Pl- 



We define the Jack symmetric polynomials whose arguments are power sums as Px'\{Pn{zj)^) = 
P\{zl+i^ ■■■,zn;P), where a = 1//?. Another important duality relation between the Jack polyno- 
mials with couplings /? and 1//? is given by 

U^a{Plr\{Pn})) = "^^Ay^^lPn}), (M) 

where c\{a) = Y\s^\iaa{s) + l{s) + 1) and c'^{a) = Yls^xicto-is) + l{s) + a). In Eq. (fT4|) . is an 
involution, an automorphism on the ring of symmetric polynomials, and is defined by 

UJaiPn) = -{-iTapn. (15) 

Using the second duality relation Eq. (I14p . we can rewrite Px' in Eq. (I12|) as 

(16) 



/3 13 

cUa) 



We should note here that the argument of P^^ in the right hand side of Eq. (jl6p is not power- 
sum Pn itself but —pn and hence -P("H{— Pn(-Zj)}) 7^ -PaC-^l+i, -2L+2, •••) -^Tv; In other words, 
^x"^ i{~Pnizj)}) is expanded by the original Jack polynomials P^{zl+i, zl+2, ■■■,zn',P) with |//| = 
|A|. By substituting Eqs. (I12p and ()16|) into Eq. ()10p . we formally obtain 

1 1 



p{wi,...,wl;zi,...,zl) = -py^o(wi, wl)^o(^i, zl) 



X E (^if({-Pn(^i)})^-Plf({-Pn(^.)}))LL ? 'fl -PA.(lZ^l,.-.,IIJL;/3)PA.(^l,.-.,^L;/?)- 



Ai,A. C^,(«)ci,(«) 

(17) 



We stress that the form of the reduced density matrix (jl7p is exact even when (N — L) 
is finite. Let us see the structure of the reduced density matrix more closely. Since 
{P'^^{{-Pn}).Ptf{{-Pn}))'N-L = ^hen |Ai| / IA2I, (dZ]) is block-diagonal and the size of each 
block is d[m) x d{m) {m = 0, 1, (i{N — L) x L), where d{m) is the number of partitions A satis- 
fying /(A) < L, 1{X') < j3{N — L), and |A| = m. Therefore, in principle, we can numerically obtain 



8 



exact eigenvalues of the density matrix by diagonalizing all the blocks in (jl7p . Although the orig- 
inal problem is reduced to the finite-dimensional eigenvalue problem, it is also difficult to evaluate 
the eigenvalues of submatrices when d{m) is large. However, if we consider (N — L) ^ oo limit, a 
considerable simplification occurs and we can evaluate the EE without any numerical calculations. 

Let us now consider the thermodynamic limit of the subsystem which traced out, i.e., (N — L) 
oo. The crucial point in our calculation is that Px^\{—Pn}) are asymptotically orthogonal with 
each other if we take the limit {N — L) ^ oo. In this limit, the reduced density matrix of our 
subsystem (jl7p becomes similar to the maximally entangled state. To see this, it is useful to expand 
the Jack symmetric functions in terms of the power sum symmetric functions [29^ as 

p 

where the power sum symmetric functions are defined for a partition p = (/0ij_£2, Pi(p)) as 

29I: 



Pp = Yli^iPpi- The coefficients Op(a) satisfy the following orthogonality relations 



cx{ay'c\{a)-'e^p{a)9^{a) = Sp^z^'a^'^^K (19) 

A 



where Zp = ni>i^'"'"^i' with mj, the number of parts of p equal to i. The coefficients Op{a) 
are nonzero if and only if |A| = From these relations, we can easily expand the power sum 
symmetric functions Pp in terms of the Jack symmetric functions as 

Pp = ^ V(^)0,^(a)c;.(a)-ip(")(K})- (20) 
p 

We stress here that the above relation itself does not depend on whether we consider the Jack 
symmetric polynomials in a finite number of variables or the Jack symmetric functions in infinitely 
many variables. By using Eqs. (jlSp and (j20p . we can formally expand Px^\{—Pn}) in terms of 
P^HiPn}) as 

Pir\{-Pn{z,)}) = CA(a)-i j;(-a)'('').,0,^(a)^,^(a)c;(a)-ip(")(K(z,)}). (21) 

p p 

Now we are ready to see the asymptotic orthoginality of P)^\{—Pn})- The scalar product of 
Px"\{~Pn}) and Px"\{~Pn}) can be represented as 

(^If ({-pn(o)}),^'lf ({-Pn(o)}))W = cxAc^r'cx.iar' 

X i-ay^''^+^^''^Zp,ZpXlia)e^plic^^^^^^^ (22) 

P1-P2 
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Suppose that the number of the particles in the subsystem traced out, {N — L), is sufficiently large, 
i.e., in the thermodynamic limit, we can simplify the scalar product in Eq. (I22p as 

lim {p(-),p(-%_^ = 5,,,,Mi(3,N-L)i^ (23) 

In this limit, we can apply the orthogonality relations ()19p to Eq. (1220 and hence we obtain 

(^'if({-Pn(0)}),<^({-Pn(0)}))W~^A,A.AA(/?,^^-L)^^^. (24) 

^Ai \^) 

We call this relation an asymptotic orthogonality of P^"\{—Pn})- The crucial point in the above 
calculation is that the sign factor (— l)'(/'i)+'(^2) -which originally comes from the expansion of 
Px"\{—Pn}) in Eq. ([22|) is canceled out. By substituting Eq. into Eq. p!7|) . the asymptotic 
form of the reduced density matrix can be expressed by the normalized basis Px as 

p{wi, ...,wl; zi, zl) ~ J2 -^a-Pa(w^i, ■■■,wl; P)Px{zi, zl; /3)^'o({w^i})*o({zi}), (25) 

A 

where Dx and Px are defined as 

D 1 -pr /?L + fl^(.)-/?r(.) 

' (^^)f-i/?^+«'(^) + i-m«) + i) ^ ^ 



and Pa = P\/ \/ {Px, P\)'l^ respectively. Then we can obtain an exact expression for the EE in the 
thermodynamic limit as 

5^,L = "Xl^AlogDA. (27) 
A 

Although this is the exact expression for the EE in the large- (A^ — L) limit, it is formidable to sum 
up all Dx log Dx because they depend on A in a complicated way. To see the physical meaning of 
this value, let us now evaluate the upper bound value of the EE. Since the reduced density matrix 
p{wi, ...,WL', Zl, zl) has already been normalized, we immediately notice that Tr p = Dx = 1. 
Under this constraint, — ^x P'x log -Da takes the maximum value when all Dxs are equal. We can 
take this maximum value as the upper bound. This maximization corresponds to neglecting the fact 
that Dx depends on the shape of the Young tableau. From the viewpoint of quantum information, 
we can say that the reduced density matrix (|25p can be approximated by a maximally entangled 
state. The upper bound value of the EE is completely determined by the number of allowed 
tableaux. Since the allowed partitions in the duality expansion Eq. ()12p satisfy /(A) < L and 
/(A') < [3{N — L), the total number of allowed tableaux is easily obtained as (^(^~^^)+^)_ Then 
the upper bound value of the EE is given by 

5.,.<5]^r = log(«^-/'+^), (28) 
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where the equahty holds when (3 = 1, i.e., the free-fermion case. Although it is one of the general 
properties that the EE is invariant under the replacement L^N — L,N — L^L, the upper bound 
itself does not satisfy this property: 5^°^'^'^ 7^ 5^°7v'i^i- This fact means that Sn,n-l approaches 
5^°^"^^, not S^°^:^^ when {N - L) 00. The upper bound 5^°^'^'^ enables us to understand 
the physical meaning of the EE in the ground state of the CS model. We now try to explain it 
in terms of exclusion statistics. In the ground state of the CS model, occupied quasi-momenta 
ki{0) are separated by /? — 1 unoccupied ones. We can schematically describe this configuration as 
Figma). In our calculation of the EE, tracing out one particle from the A^-particle ground state 
corresponds to the decimation of one quasi-momentum from the Fermi sea. In other words, one 
1 is removed from the Fermi sea when we trace out one of coordinates Zj. As we said before, /3 
quasiholes (/3 zeroes) are created in the Fermi sea in this process (see Fig. [U^b)). It is now obvious 
that tracing out (A^ — L) particles from the ground state corresponds to the decimation of (A^ — L) 
quasi-momenta from the Fermi sea and the creation of P{N — L) quasiholes in the Fermi sea. The 
number of possible intermediate states consisting of L particles and P{N — L) quasiholes can be 
counted as follows. First we recall that the Fermi sea consists of A^ I's and (/? — 1)A^ O's. After 
the decimation of the (A^ — L) quasi-momenta, the configuration of the state consists of L I's and 
(/? — 1)A^ + (A^ — L) O's with the exclusion constraint such that any two I's are separated by more 
than (P — 1) O's. Finally, we notice that the number of possible intermediate states is identical 
to that of possible configurations of I's and O's with the constraint and can be easily obtained as 
• Here we can see that the upper bound of the EE 5^°)^'^'^ is equal to the logarithm of 
this number. It is also remarkable that 5^°)^'^*^ coincides with the upper bound value of the EE in 
the Laughlin state if we identify m = /?, where m denotes the inverse of the filling factor |26l |. 



It would also be possible to interpret S^'^'^ in terms of the flux attachment in the context of the 
quantum Hall effect. While 5^°)^'^*^ provides a natural way to understand the EE in the CS model 
in terms of the fractional exclusion statistics, we can also obtain a more accurate value of the EE 
by taking it into account that Dx depends on the shape of the Young tableau A. Comparing this 
value with S^f^^'^, we notice that the subleading term, Sn,l — •S^f^^'^, does not depend on the 
total number of particles A^ but only on the coupling /3 and L. A similar universal property has 
already been found in the study of the one-particle EE of hard-core anyons on a ring, where the 
subleading term depends only on the anyonic parameter 9 [30]. The details of the calculations and 
the difference between Sn^l and 5^°)^°'^ in the thermodynamic limit are argued in Appendix El 
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IV. SUMMARY AND DISCUSSIONS 



In this paper, we have studied the entanglement entropy between two blocks of particles in 
the ground state of the Calogero-Sutherland model. We have obtained the exact expressions for 
both the reduced density matrix of the subsystem and entanglement entropy in the limit of an 
infinite number of particles traced out. In our calculation, the duality relations between the Jack 
symmetric polynomials with coupling j3 and those with 1//3 have played a crucial role. From the 
obtained results, we have estimated the upper bound value of the EE by a variational argument. 
We have also found that the upper bound value itself has a clear physical meaning in terms of 
fractional exclusion statistics. This interpretation indicates that entanglement between subsets 
of particles enables us to extract interesting properties in a wide range of systems with fractional 
exclusion statistics. It is also remarkable that this upper bound coincides with that of the Laughlin 
state in fractional quantum Hall systems when we identify the inverse of the filling factor m = (5. 

While we have studied the EE between two blocks of particles, it would also, of course, be 
important to study the EE between two spatial regions in the ground state of the CS model. In spin 
systems on a lattice such as the XY spin chain in a transverse magnetic field, it is possible to perform 
an exact analysis of the EE between two spatial blocks with the aid of the Fredholm determinant 
technique Q]. This technique based on the Riemann-Hilbert problem also plays a crucial role 
in the computation of the correlation functions for random matrices. On the other hand, it is 
known that the CS model is identical to Dyson's brownian motion model of the circular ensembles 
with /3 = 1,2,4 



31[. Thus, it is promising to obtain the EE for spatial partitioning by applying 
the Fredholm determinant technique. It would also be interesting to investigate entanglement 
properties injntegrable lattice models with inverse square interactions such as the Haldane-Shastry 
model 



33( 1 and the long-range supersymmetric t-J model [34,]. It remains an interesting issue 



whether our method developed in this article can be directly applied to these systems by using the 



freezing trick 



351,1331. 
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APPENDIX A: MORE DETAILED ANALYSIS 

In this appendix, we discuss the more detailed analysis of the EE (j27p and the universal sub- 
leading correction of the EE. Although both Sn^l and 5^°)^'^'^ go to infinity in the thermodynamic 
limit: N ^ oo and L is fixed, the difference S^'^'^ — Sn^l is finite. The strategy to show this 
fact is to rewrite the sum over partitions as the integral over continuous variables. This method 
is similar to the calculation of the dynamical correlation functions in the CS model by Lesage, 
Pasquier, and Serban [21]. Let us start with rewriting Eq. ([26]) in terms of parts of A: 

1 J^ T{f3{L-j) + l)T{Xj + P{L-j + l)) _ P^Ll{P{N-L)y. J^ T{Xj + P{L-j + l)) 

' © M - j + i))r(Ai + PiL - j) + 1) my- /i r(A,- + piL - j) + 1) 

(Al) 

Introducing the new scaled variables tj = Xj/N and using the Staring formula: r(x + 1) 



V27rx(x/e)^, we obtain the simple expression for Dx- 



—f{t^,...,tL;l3), (A2) 



where /(ti, . . . ,tL', P) = Y[f=i N —>■ oo, we can replace the sum over {Xj} with the integral 

over {tj}: 

^ E [jt,...dtL, (A3) 

0<A£,<-<Ai</3(Af-L) 



D 



where D is the region satisfying < t/, < • • • < ti < /?. From these results, Tr p = '^x^>^ 
evaluated as 

L\ f 1 



1 



/3 



dtt^^-^] =1. (A4) 



L 



This is consistent with the normalization condition of p. Similarly, the EE Sn^l can be rewritten 
in terms of the integral over {tj}: 

= L log iV - log L! + (/3 - 1)L log /? - " " " '^^'"^ ^' ^^^^ 
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We can exactly evaluate the integral of the last term as follows: 

r-/3 rP 



f dti--- dtLf log / = i- / dh-- - [ dtLih ■ ■ ■ tLf-^ V log t^r^ 

JD ^- Jo Jo ~[ 

(/5 log /3 - log + /?-')■ (A6) 



(L-1)! 

Thus Sn,l = LlogN- logL! + L(l - P'^). On the other hand, since = log ^ 

L log N — log L! + L log P, we finally obtain 

S^T"" - ^N,L ~ ^(log P-l + p-'). (A7) 

Therefore, the subleading term of the EE does not depend on the total number of particles but 
only on L and the coupling of the CS model /?. Note that the right hand side of Eq. ([AT]) vanishes 
only for /? = 1. This result means the EE can saturate the upper-bound entropy only for the free 
fermion case. 
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